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Abstract 

The statistical properties of velocity gradients in a wall-bounded turbulent channel 
flow are discussed on the basis of three-dimensional direct numerical simulations. 
Our analysis is concentrated on the trend of the statistical properties of the local 
enstrophy uj'^{x,t) and the energy dissipation rate e{x,t) with increasing distance 
from the wall. We detect a sensitive dependence of the largest amplitudes of both 
fields (which correspond with the tail of the distribution) on the spectral resolution. 
The probability density functions of each single field as well as their joint distribution 
vary significantly with increasing distance from the wall. The largest fluctuations 
of the velocity gradients are found in the logarithmic layer. This is in agreement 
with recent experiments which observe a bursting of hairpin vortex packets into the 
logarithmic region. 
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1 Introduction 



Turbulence is associated with large fluctuations of velocity gradients which 
appear preferentially at the smallest scales of the flow. The amplitudes of the 
fluctuations exceed the mean values by orders of magnitude when the Reynolds 
number of the flow is sufficiently large. This behaviour is known as small-scale 
intermittency p^. It has been discussed recently that a deeper understanding of 
fluid turbulence as a whole requires a detailed resolution of the intermittent 
dynamics at the small-scale end of the inertial range [2]. This will include 
scales that are smaller than the (mean) Kolmogorov dissipation length tjk = 
v^l^j^e)^!^ with the kinematic viscosity v and the mean energy dissipation 
rate (e). Although significant progress in measurement techniques has been 
made, the finest structures remain still spatially unresolved |3] or the flows are 
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limited to very low Reynolds numbers [1]. Direct numerical simulations can 
reach today sufficiently high Reynolds numbers while resolving the gradients 
at the small-scale end of the inertial range cascade properly [51IS1I7] . 

The small-scale structure and statistics of turbulence has been mostly studied 
for the case of homogeneous, isotropic and statistically stationary turbulence. 
A wall-bounded shear flow consists of a boundary layer which is dominated by 
coherent streamwise structures and a central bulk region in which they are ba- 
sically absent. The statistics in the wall- normal direction is inhomogeneous and 
requires a height-dependent analysis on account of the varying strength of the 
mean shear. As a consequence, there is significantly less work on this aspect for 
wall-bounded flows than for homogeneous isotropic turbulence [8]|9]|10llll] . The 
coupling between large-scale shear and intermittency in the inertial cascade 
range has been studied in channel flow simulations [T^HS]- Further studies in 
the same flow configuration demonstrated the violation Kolmogorov's refined 
similarity hypothesis (K62) [14j when approaching the walls [15|IT6] . It was 
therefore suggested to replace the third-order increment moment in K62 by 
a second-order one which is directly coupled to the large-scale shear in the 
corresponding balance equations. As a consequence, shear-improved eddy vis- 
cosities in Smagorinsky subgrid-scale models can improve the representation 
of Reynolds stresses as shown in Ref . . 

In comparison with homogeneous, isotropic and statistically stationary turbu- 
lence, the presence of an inhomogeneous direction increases the level of flow 
complexity significantly. With regard to the velocity gradient statistics one can 
identify the following two aspects. First, it is more challenging to measure all 
nine derivatives of the velocity gradient tensor with the sufficient resolution 
in an open shear- flow-setup (see e.g. [IB])- Second, it is frequently believed 
that shear flow turbulence is in a state of local isotropy at the small-scale 
end for larger Reynolds numbers [9] . More recent experimental and numerical 
studies demonstrated however that significant deviations persist, in particular 
when higher-order moments are discussed p^20f21] . It can thus be expected 
that some fundamental differences to the homogeneous isotropic case persist 
to very high Reynolds numbers. All this suggests to our view a systematic 
study of the height-dependence of the statistics of the velocity gradient fields 
in turbulent shear flows. 

In the present work, we want to make a ffist step in this direction and con- 
duct an analysis of the small-scale statistics of the velocity gradient in a wall- 
bounded shear flow. It is based on three-dimensional direct numerical simula- 
tions of a turbulent channel flow. Our study will be focussed to two fields, the 
local enstrophy field and the energy dissipation rate field. The gradient tensor 
of the velocity fluctuations, rriij, can be decomposed into a symmetric and an 
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Fig. 1. Sketch of the channel flow geometry and the coordinates. The simulation 
domain is x L„ x Lj; = 47r x 27r x 2 in units of the half-width L. 



antisymmetric part 



(1) 



dxj 



The local enstrophy field, uj"^ = (V x v')^, probes the magnitude of the an- 
tisymmetric part ttij = —eijkUJkf^ of the velocity fluctuation gradient tensor. 
Here eijk is the fully antisymmetric Levi-Civita tensor. The energy dissipation 
rate field, e = 2vsijSji, measures basically the magnitude of the symmetric 
part Sij of iTLij. The ensemble average of both fields are connected by 



for homogeneous isotropic turbulence. This does however not imply that the 
local statistical fluctuations are synchronized [3f22] . In addition to homoge- 
neous isotropic turbulence, we have to conduct the statistical analysis for both 
fields as a function of distance from the wall and to separate the turbulent 
fluctuations from the mean flow. 

The outline of the paper is as follows. We discuss the equations of motion 
and the numerical model in the next section. Afterwards we investigate the 
dependence of the small-scale statistics on the spectral resolution which sets 
the stage for the analysis. The subsequent section investigates the single-field 
and joint statistics of the local enstrophy and dissipation rate fields as a func- 
tion of the distance from the wall. We summarize our findings and give a brief 
outlook at the end. 




(2) 



2 Equations and numerical model 



The numerical method solves the Navier-Stokes equations for an incompress- 
ible flow which are given in dimensionless form by 

^ _ X = - V (p + t;V2) + — V'-y, V-v = 0, (3) 

with ri = V X t; and boundary conditions 

v = Q at z = ±1. (4) 

Figure [T] illustrates the geometry and the coordinates of the channel flow. The 
applied pseudospectral method is based on Fourier series in the horizontal 
directions x and y and a Chebyshev polynomial expansion in the vertical z- 
direction [23112^ . The solenoidal velocity fleld is represented in terms of two 
scalar quantities if) and using the poloidal-toroidal decomposition 

v{x,t) = V X (V X ezip{x,t)) + V X ezip{x,t). (5) 



Equations for (p and ifj are derived by taking the curl and twice the curl of 
the momentum equation and projection onto the vertical direction. We obtain 
two equations for the vertical velocity component Vz = —^hf and the vertical 
vorticity component Qz = —^h'^P, where = + dy denotes the horizontal 
Laplace operator. The quantities Vz and Qz determine the velocity fleld up 
to a mean flow U{z,t)ex + V{z,t)ey. Equations for U and V are obtained 
by averaging the momentum equation over horizontal cross-sections of the 
periodicity domain. Upon introducing the deflnitions 

C = ^lz, ri = V^t;^, i = Vz (6) 



for ease of notation, the evolution equations based on the poloidal-toroidal 
representation take the form 



V^C - Re^tC = -Ree, ■ [V x x Q,)] , 

V^r] - Re dtT] = Re [dzV ■ {v x fl) - ez ■ {v x fl) 

v'e = V, 

dlU - Re dtU = Re [dz (v^Vz)^ + d^p] , 
djV-RedtV = Re . 



(7) 

(8) 

(9) 
(10) 



The angular brackets {■)a denote horizontal averages. The mean flow equa- 
tions also contain the mean driving pressure gradient d^^p- The appropriate 
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Fig. 2. Vertical profiles of the Reynolds shear stress ris (left) and the energy dissipa- 
tion rate e (right) averaged in planes of constant z. The wall-normal coordinate z is 
given in wall units. The results of small- (128^ points, blue curves), medium- (256^ 
points, green curves) and high-resolution (512^ points, red curves) simulations are 
shown. In each case, the quantities are made nondimensional with the wall friction 
velocity Ur from the simulation with 512^ points. 



boundary conditions are readily derived using the incompressibility constraint 
V ■ i; = 0. We obtain 



^ = ( = 5^^ = t/ = 1/ = at z = ±l. 



(12) 



For the time discretization of the evolution equations we use a method of 
second-order accuracy. If we write such an equation symbolically as dtf = 
Lf + .Afif) where C denotes a linear operator and A/" the remaining terms, 
our time-stepping scheme reads 



2At 



(13) 



where At is the time step. The left hand side approximates the time derivative 
dtf at the time level n + 1 using the previous two time levels. The linear term 
Cf is treated implicitely, and J\f explicitly through the second-order Adams- 
Bashforth method, where the prefactors correspond to a linear extrapolation 
to the time level n+1 [25] . The solution procedure for the coupled eqns. ([8]) 
and ([9]) requires auxiliary functions to satisfy the boundary conditions d^^ = 
for eq. ([8]) as explained in j25j. The numerical code implements distributed- 
memory parallelization through MPI. 
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3 Single quantity statistics 



3.1 Spectral resolution effects 

Our statistical analysis is focussed to the small-scale properties of the channel 
flow turbulence. In order to set the stage for this investigation, we flrst compare 
statistical properties at three different spectral resolutions and the same set of 
parameters. The flrst case with 128^ (or more precisely 128'^ x 129) grid points 
corresponds with the classical simulatioij^ of Kim, Moin and Moser [26] in 
terms of the vertical resolution. The second case is two times better resolved 
in each space direction with 256^ grid points. The highest spectral resolution 
is four times better resolved than the flrst one and contains 512^ (or more 
precisely 512^ x 513) collocation points. The three runs were performed with 
identical flxed volume flux such that the resulting friction Reynolds number 
approaches the grid-independent limit Rcr = 180. The deflnition 

„ vLtL . , 

Rer = ^. (14) 
f 

is based on the wall friction velocity m,-, the half-channel height L, and the 
kinematic viscosity by v. We follow the usual convention and denote the tur- 
bulent velocity fluctuations by 

v'{x,t) = v{x,t) -U{z)e^. (15) 

Figured (left) compares the Reynolds shear stress ri3 = {v'^v'^) function 
of z'^ = zv/ur- The three vertical proflles of ris coincide quite well, except for 
slight deviations around the maximum in case of the lowest resolution. The 
additional subscript t in {■)A,t indicates time averaging. We repeat the analysis 
for the vertical proflle of the energy dissipation rate averaged over time and 
in planes at constant distance from the wall. The quantity is deflned as 

1/ I ( dv' dv'A\ 
ea„{z) = {eix,t))A,t = 2iy{s,,s,,)A,t = 2 ^ 8^) / ' ^^^^ 

The plot yields almost perfectly collapsing curves for most of the channel. De- 
viations for the smallest resolution exist for the viscous layer at z^ < 20. We 
can conclude that the lowest resolution is sufficient when flrst-order moments 

^ More precisely, the grid resolution in [26] was x Ny x Nz = 190 x 129 x 160 
(with y being the wall-normal coordinate) for a channel with the same sidelengths 
as in the present work. 
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Fig. 3. Comparison of the probability density functions of the energy dissipation rate 
(upper row) and local enstrophy (lower row) for the two spectral resolutions. Data 
are obtained in planes at four different distances from the wall which are indicated 
in the legends. The left column stands for the lowest resolution case, the right one 
for the highest resolution. Both quantities are given in units of their statistical mean 
in the corresponding plane. We denote Cavr = {^)A,t and w^^, 



are considered, such as shear stresses. It is therefore suitable for many pur- 
poses, such as validations of closures of the Reynolds averaged equations [27] . 
We can however see already, that the mean dissipation rate profile deviates 
significantly in the viscous layer. 

A full statistical information requires the knowledge of all moments or the 
probability density function (PDF). We proceed therefore by comparing the 
PDFs of the energy dissipation rate, e{x, t), and the local enstrophy, uj'^{x, t) = 
(V X f')^, for both resolutions. Four different distances have been chosen: 
= 10 is in the viscous buffer layer, z+ = 30 at the begining of the loga- 
rithmic layer where the Reynolds shear stress has its maximum (see Fig. [2] 
(left)), = 72 in the logarithmic region, and = 180 in the midplane 
of the channel. While the first-order moments are basically unchanged at the 
higher resolution, we can detect clear differences in the far tails of the PDFs for 
all cases. The far tails are exactly the parts of the PDF that become relevant 
for higher-order moments. Their sufficient resolution requires the efforts made 
in the present work (see also [S] for a detailed discussion). 
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Table 1 

Moments of the normalized energy dissipation rate as a function of height and 
resolution. 





<e2> 


< > 


< > / < e2 >2 


z+ 


128 256 512 


128 256 512 


128 256 512 


10 
30 
72 
180 


2.12 2.15 2.14 
2.25 2.93 3.14 
2.46 2.96 3.15 
2.55 2.81 3.06 


7.73 7.93 7.83 
9.92 17.95 21.83 
13.00 20.43 24.36 
14.56 18.67 23.07 


9.59 9.45 9.34 
15.45 20.40 25.44 
21.34 32.53 36.77 
25.24 30.47 35.94 



3.2 Height dependence of statistics 

Let us turn now to the high-resolution case of both fields in the right column of 
Fig. [3l Note first that the amplitudes are always rescaled by the corresponding 
mean values for the given plane at fixed z'^ . Close to the wall, the turbulent 
velocity fiuctuations are small. The mean dissipation profile indicates however 
that this is not the case for the magnitude of the derivatives. As we can see 
in Fig. [2], the mean energy dissipation rate is largest close to the wall and 
drops by more than two orders of magnitude toward the channel centre plane. 
The tails of both distributions are the sparsest ones in the vicinity of the 
wall. While the means are large the variations about the means remain small. 
The extension of the tails varies with increasing distance from the wall as the 
data for the three other values of indicate. Since the mean values decrease 
monotonically towards the midplane, high-amplitude events are shifted ever 
further into the tail. The fattest tails are observed at z^ = 72, i.e. in the middle 
of the logarithmic region. Particularly here, high amplitudes can be expected 
further away from the wall due to bursting events. They are for example caused 
by coherent hairpin vortices that can grow out into the logarithmic layer. 
Such a large-scale hairpin vortex packet formation has been recently suggested 
by Adrian and co-workers on the basis of experiments and direct numerical 
simulations pjj. Our statistical analysis is consistent with this observation. 

The high resolution case reproduces another feature which is well-known from 
the homogeneous isotropic case. The tail of the local enstrophy is fatter than 
that of the energy dissipation rate [3|6] . Simulations of isotropic box turbulence 
by Donzis et al. [6j suggest that this difference in the tails decreases with 
increasing Reynolds number, a point that certainly has to be studied in the 
shear fiow well. 

Height and resolution dependences of higher moments of the dissipation rate 
are given in table [H The fatter tails of the PDFs at higher resolution account 
for the growth in the second and third moment with the resolution for all z^ 
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values except in the viscous region at = 10. Likewise, the flatness (e^)/ (e^)^ 
is growing with the resolution. This quantity is also indicative of the degree 
of intermittency at the smallest lengthscales. Remarkably, it is not constant, 
but increases significantly from = 10 to = 72, i.e. the level of small- 
scale intermittency depends on the wall distance outside the viscous region. 
However, it changes only slightly from = 72 to the midplane. 

Regarding the relation between planar averages of enstrophy and energy dis- 
sipation we remark that they are connected by the relation 

{e)A = iy{u^')A + iydl{v'J')A (17) 



in a wall-bounded shear fiow. In contrast to Eq. ([2D, this relation contains 
an additional derivative term arising from the averaging which is restricted to 
planes of fixed z. We found that the second term on the right hand side of f|T71) 
is at least two orders of magnitude smaller than {u)'^)a- We have also verified 
that relation f[T7l) is satisfied in our numerical simulations. Furthermore, the 
statistics of the energy dissipation rate field deviated from the log-normal 
distribution for all four planes, particularly in the far tails at the small and 
large amplitudes. This property is similar to isotropic turbulence. 



4 Joint statistics of energy dissipation and local enstrophy 



In view of the differences in the PDFs for enstrophy and energy dissipation 
it is interesting to study how both fields are statistically correlated. For this 
reason we plot in Fig. IHthe joint PDF of e and uj"^ normalized by both single 
quantity PDFs 

n(e...^) = ^!f^ (18) 

p{t)p[u^) 



High amplitudes of H indicate a strong statistical correlation between both 
fields. The PDF of the data in the midplane has the broadest support. Very 
high-amplitude events of dissipation and vorticity are strongly correlated as 
indicated by the local maximum which is stretched out into the upper right 
corner of the plane. This maximum is even more pronounced in the logarithmic 
layer as can be seen in the upper right panel of Fig. HI The local maximum of H 
in the outer right corner becomes narrower with decreasing distance from the 
wall. It is in line with an overall decrease of the support for highly correlated 
events. The observed trend is consistent with our observations for the single 
quantity statistics. Again, we can detect a decrease of the fiuctations of both 
fields about the mean. As a guide to the eye, we added the line e = vu)"^ 
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Fig. 4. Joint and normalized PDFs for the energy dissipation rate e and the square of 
vorticity as given by Eq. ([18]). Data are taken again at four z+-positions. Upper 
left: = 180. Upper right: z"*" = 72. Lower left: z'^ = 30. Lower right: z'^ = 10. 
Results of DNS with 512^ grid points are shown, the joint PDF fields are normalized 
by the product of separate PDFs for ejeavr and uP' jio^^^. The white diagonal line 
in the lower right figure corresponds with e = vijp' . The colour coding is the same 
in all four figures and given in units of decadic logarithm. It varies from 10° (red) 
to 10""^ (blue). 
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Fig. 5. Isosurfaces of the energy dissipation rate (upper figure) and the local enstro- 
phy (lower figure). The isolevels are taken at three times the mean dissipation and 
the mean squared vorticity. 

(in white) to the lower left panel, i.e. to the joint-PDF plot closest to the 
wall. It can be seen that a broad range of amplitudes of both fields is basically 
concentrated now around this line, which connects the means of both fields (see 
eqns. ([2]) or ( |T7|) ). To conclude, at all distances from the wall there is a strong 
statistical correlation between the local enstrophy and the dissipation. Close 
to the wall high-amplitude events of both fields seem to become synchronized. 
The support of the statistically correlated events varies significantly across the 
channel. 

Figure [5] displays the structures of both fields. The largest amplitudes of both 
fields are clearly concentrated close to the walls. We also observe the hairpin- 
like vortices that burst into the channel center (see lower panel). Isosurfaces 
of the energy dissipation rate display more sheet-like structures (see upper 
panel). They reach also out into the logarithmic region. Such shear layers 
can be formed when pairs of streamwise streaks are lifted up |28] or when 
they become linearly unstable. The inspection of the structures in both plots 
supports thus our statistical results obtained before. Maxima of both fields 
are found close together in space. In particular, we expect, that exactly the 
bursting events will contribute significantly to the local maximum of the joint 
PDF in the upper right panel of Fig. |H A firm proof would require a separation 
of these coherent structures in the data record, which will be part of our future 
work on this subject. 
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5 Summary and conclusions 



We have presented a study of the statistics of the energy dissipation rate and 
local enstrophy in a turbulent channel flow. We showed that a full statisti- 
cal analysis of both fields requires a spectral resolution which is larger than 
the standard one. This was one objective of our study, namely to make this 
case in a turbulent wall-bounded shear flow. It clearly limits the range of 
accessible Reynolds numbers. Based on this result, we investigated the height- 
dependence of the statistics of both fields in a first example. We showed that 
the variations of both fields about their mean values are strongest in the 
logarithmic layer. They decrease slightly towards the channel centre and sig- 
nificantly towards the wall. The statistics of the velocity gradients is thus 
strongly height-dependent. It is clear that a more comprehensive parametric 
study has to follow, e.g., in order to analyze these trends at different Reynolds 
numbers. 

There is a further interesting aspect onto which our studies can shed new 
light. Recently, new ideas for multifractal subgrid-scale models have been de- 
veloped which are based on exactly such statistical fluctuations of velocity 
derivatives which have been studied here [21]. The model of has been 
successfully validated in homogeneous isotropic turbulence and it seems to be 
a hopeful candidate for application to shear flow turbulence. Investigations in 
this direction are under way. 
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